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SUPPLEMENT TO: A COMPUTER PROGRAM FOR THE DESIGN AND ANALYSIS OF 

LOW-SPEED AIRFOILS 


Richard Eppler* and Dan M. Somers 
Langley Research Center 


SUMMARY 

Three new options have been incorporated into an existing computer 
program for the design and analysis of low-speed airfoils. . These options 
permit the analysis of airfoils having variable chord (variable geometry), a 
boundary-1 ayer displacement iteration, and the analysis of the effect of 
sinqle roughness elements. All three options are described in detail and are 
included in the FORTRAN IV computer program which is available through COSMIC 


INTRODUCTION 

A conformal -mapping method for the design of airfoils with prescribed 
velocity-distribution characteristics, a panel method for the analysis of the 
potential flow about given airfoils, and a boundary- layer method have been 
combined. With this combined method, airfotls with prescribed boundary-layer 
characteristics can be designed and airfoils with prescribed shapes can be 
analyzed. All three methods and the FORTRAN IV computer program for the 
numerical evaluation of these methods are described in reference 1. 

Three new options have been incorporated into the computer program 
described in reference 1. The previous version of the program (ref. 1) was 
capable of analyzing an airfoil with a simple flap. In the present version, 
an option has been added which allows the analysis of an airfoil having 
variable chord (variable geometry). The method of reference 1 did not contai 
a boundary-layer displacement Iteration. An iteration procedure has been 
included in the present version. The third option to be added permits the 
analysis of the effect of single roughness elements. The input for all three 
options is described in detail. 

Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


♦Professor, University of Stuttgart, Stuttgart, West Germany. 



SYMBOLS 


Values are given in SI units. 

Cf boundary- layer skin-friction coefficient 

c airfoil chord, m 

cy section profile-drag coefficient 

c^ section lift coefficient 

section pitching-moment coefficient about quarter-chord point 
h height of roughness element normal to surface, m 

Is lower surface 

R Reynolds number based on free-stream conditions and airfoil chord 

R. Reynolds number based on local conditions and height of roughness 

n element 

U potential -flow velocity, m/s 

U„. free-stream velocity, m/s 

u. x-component of velocity In turbulent boundary- layer at height of 

n roughness element, m/s 

us upper surface 

v local velocity on airfoil, m/s 

x airfoil abscissa, m; axis in streamwlse direction, tangential to 

surface 

x R chord location of roughness element, m 

y airfoil ordinate, m 

a angle of attack relative to zero-lift line, deg 

A incremental change in quantity 

6 1 boundary-layer displacement thickness, m 

6 2 boundary-layer momentum thickness, m 

v kinematic viscosity, m^/s 
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p 


air density, kg/im* 

2 

t q shear stress at wall, kg/m‘S fc 

PROGRAM AVAILABILITY 

The program Is available at a nominal fee through the following 
organization: 

Computer Software Management Information Center (COSMIC) 
112 Barrow Hall, University of Georgia 
Athens, Georgia 30602 

Request the program by the designation PROFILE LAR-12727. 


VARIABLE GEOMETRY 

The previous version of the computer program (ref. 1) allowed the shape of 
an airfoil analyzed by the panel method to be altered so as to correspond to 
the deflection of a simple flap. Thus, that version only permitted the rota- 
tion of a portion of the airfoil, the flap, about a specified hinge point. 
Chord-increasing flaps were not allowed. The present version of the program 
can analyze this form of variable geometry. It should be noted that, while the 
airfoil shape which results from the exercise of this option does have an 
Increased chord, it does not contain a slot and, thus. Is still a single- 
element as opposed to a multi-element airfoil. An application of this capa- 
bility Is described in reference 2. 


FLAP Card 

The variable-geometry option Is selected by setting NUPU * 1, 2, 3, or 4 
on the FLAP card. 

NUPA, NUPE, and NUPI are neglected. 

NUPU = 1 — The F-words specify the points to be deleted. The five digits 
of F.j are denoted aaabb. Points aaa through aaa + bb 
are deleted. If bb * 00, only point aaa is deleted. 

It Is recommended that the F-words be specified with decreasing values of 
aaa as the points after point aaa (higher point number) are renumbered. This 
means that aaa for F, should be greater than aaa for F 2 which should be 
greater than aaa for Fj and so on. 

Only one FLAP card with NUPU =1 Is allowed. 
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NUPU * 2 - The F-words specify points to be added to the upper surface. 

The new points are added after the point on the upper surface 
having the greatest x/c remaining after the deletions which 
resulted from the FLAP card with NUPU * 1. Thus, if point 1 
(x/c ■ 1) was not deleted, only points with x/c > 1 can be 
added. 

0.01F 1 = Xj/c [f5.4J 

0.01F 2 * y|/c [f5 . 4 ] 

(F 3 .F 4 ) = (x 2 /c,y 2 /c) and so on 

It should be noted that the new points must be In order of Increasing x/c. 

NUPU = 3 — The F-words specify points to be added to the lower surface. 

The F-words are interpreted just as they are for a FLAP card 
with NUPU * 2. 

NUPU * 4 - The F-words specify additional points to be spllned In between 
the points available so far. The F-words are interpreted just 
as they are for an FXPR card. (See ref. 1, p. 45.) 

It should be remembered that the points are renumbered during the 
execution of each of the preceding FLAP cards. 

The panel method Is called automatically after a FLAP card with NUPU = 4 
Is read. Following this card, any other cards (in the- proper sequence, of 
course) are allowed except another FLAP card. Only airfoil coordinates gener- 
ated in the design mode or read In following an FXPR card can be altered by 
FLAP cards with NUPU = 1, 2, 3, and 4. Thus, a FLAP card with NUPU =1,2, 

3, or 4 cannot follow another FLAP card. 


Example 

The following card sequence illustrates the use of the variable-geometry 
option. 
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TP*1 


Af 4 i;?«50 ^00 1450 000 1650 9 00 1950 400 9050 500 ^950 600 2450 700 

6*4 ?*50 900 9940 900 9090 1090 000 1*00 9950 100 9550 900 4950 900 

min » u m * i hTb maa uTTh m>mnnmmMiiiMnin»ii>iMiHBii»iiiinMM awn n n aj 

AA4 4550 900 4750 100 4950 000 5150 9500 6000 400 


T*9! 


»«t ij 


T9A1 


M « « 


TR99 


mm muimN 


664 400 1950 900 1000 7?0 400 950 200 900 710 100 90Q 000 000 


rTTSji '» m»n»»»iiBiiiMiiii»>i»»M»BnmitiiMHM>iii»MHimin>MniiMii ■»» i»n» » ■] 


OLT6 


1 1 000 


nBBBn«awBBHBBB*i»Baa im aBwaqBaaiiBBBMBaBBBWBBBBBBBiBpaBBBaaaB 


ilia n a a a a a 


01 *G 


ilia n n i»h a a ti a a 


ansa a a a w a a a a a a a a a » a a m m a a a miimiii »»m a a h a a a a au urn a a a it a a n a a n a ijj 


FLAP 


1 5909 700 500 901 


rrrjjTWa a~i a a Vi I gawarBBBwaaaaaiiaaaiiaaaa ana iauniitMHM«HBBi»innanimiMnai») 

911000 -95019000 -900 


FLAP 


i i i an unaa auaiTaii aaar. anaaaa BaMaB»aaaiiaawaiia^awaaaaBBaawB«aaawa an n wa aan aTg 


FL^P 


3 Q400 -*65 9600 -9901 0900 -49019000 -900 


i * i 

FL9° 

r '“r 


i i t an nnwnaiMiwaiinaaaariaaawaaaaMiiaaaiiaBaaa«iaaB>iBBwaanBaaaaaaaaaaaawaawaaHa Fa] 


45940051 9 005 090 0492 000990001 900 


PLFP 


i i i a n n n a n » » a n a n o a a a a n a a a n » 

1 1 000 


i i i a it ti n a n a n a a aiiaa nannm 


aaaaauaaaiigaaaaiiaaaiiM a a a ana«BaBuaaaaaaawaaaaanaaa| 


019^5 


a a a a n a a a *i a o a a a « a a a n a a a a a n a a a a a o a a a i> a a a h it a » a a n a a a) 


i i i a ii it ii m is a ii a a a naaaaanaaawa 


rMT’F 


n i a ii n n w a a n’a a a n a a a a aa <i a a a n a a a a a n a a a a a u a a a m a a n n a a a a n a a aj 


The first FLAP card deletes points 52 through 61 as well as points 7, 5, 
3, and 2 (in the x-y-v listing, N * 51 through 60, 6, 4, 2, and 1). If the 
chord is to be increased, some of the points near the trailing edge should be 
deleted. In other words, a short distance between points is required near the 
new trailing edge, not the old one. The second FLAP card specifies two points 
for the extension of the upper surface: (x/c * 1.1000,y/c = -0.0350) and 
(x/c = 1.2000,y/c * -0.0900). The third FLAP card specifies four points for 
the extension of the lower surface: (x/c * 0.8400,y/c * -0.0265), 

(x/c « 0.9600,y/c * -0.0290), (x/c * 1.0800,y/c * -0.0490), and 
(x/c * 1.2000,y/c = -0.0900). The fourth FLAP card Inserts In the equiangular 
spacing mode (ref. 1) four points between points 52 and 53, two points between 
points 51 and 52, two points between points 50 and 51, two points between 
points 49 and 50, two points between points 2 and 3, and three points between 
points 1 and 2. The panel, method is called automatically after the fourth 
FLAP card. 


This card sequence plot into one 
tions for airfoil 664 with and without 
Two velocity distributions, each at o 
plotted. The following x-y-v listings 


diagram (fig. 1) the velocity distribu- 
te variable-geometry flap extension. 

■ 0° relative to the chord line, are 
are also generated. 
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AIRFOIL 664 16*631 0*00 

NX Y VELOCITY 01 $TRI 3UTI QNS FOR THE ABOVE ANGLES OF ATTACK RELATIVE TO THE CHORO LINE 


0 

1.00000 

0.00000 

.773 

1 

.99623 

* .00092 

.799 

2 

.98557 

.00391 

.862 

3 

.96923 

.00881 

.934 

4 

.94774 

.01491 

.970 

5 

.92110 

.02193 

.997 

6 

.88964 

.03005 

1.020 

7 

.85407 

.03927 

1.048 

8 

.81512 

.0494? 

1.001 

9 

.77353 

.06020 

1 .121 

10 

.73C08 

.07122 

1.167 

11 

. 6e 549 

.00197 

1.220 

12 

.64043 

.09167 

1.281 

13 

.59497 

.09937 

1.318 

14 

.54869 

.10402 

1.326 

1? 

.50167 

. 1 OS 40 

1.331 

16 

.45437 

.11029 

1.334 

17 

.40727 

.11060 

1.335 

18 

.36097 

.10939 

1.335 

19 

.31564 

.10670 

1.334 

20 

.27205 

•10262 

1*330 

21 

.23051 

.09720 

1.324 

22 

. 19145 

.09055 

1.315 

23 

.15521 

*08?77 

1.301 

24 

.12216 

.07401 

1.281 

25 

.09258 

.06441 

1.252 

26 

.06674 

.05416 

1.211 

27 

.04407 

•04348 

1.150 

28 

.02714 

.03261 

1.056 

29 

.01371 

.02183 

.909 

30 

.00468 

.01155 

• 661 

31 

.00023 

.00229 

.197 

32 

•00146 

-•00521 

.562 

33 

.00903 

-.01173 

.842 

34 

•02234 

-.01817 

.954 

35 

.04097 

-.02423 

1.004 

36 

.06471 

-.02979 

1.028 

37 

.09334 

-.0348? 

1.039 

38 

. 12641 

-.03936 

1.046 

39 

.16380 

-.04341 

1.052 

40 

.20474 

-.04693 

1.057 

41 

.24882 

-.04990 

1.060 

42 

.29552 

-.05229 

1.064 

43 

.34429 

-.05406 

1.066 
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AIRFOIL 66 A 16.63T 0*00 

H X Y VELPCITr DISTRIBUTIONS FOR 


44 

.39452 

-.05522 

1.070 

45 

•44556 

-'.05572 

1.074 

46 

.49678 

-.05546 

1.079 

47 

.54754 

-.05433 

i.oe4 

48 

•59719 

-.05219 

1,090 

49 

.64512 

-.04967 

1.096 

50 

•69117 

-•04322 

1.076 

51 

.73561 

-.03567 

1.036 

52 

.77909 

-.02623 

.957 

53 

.82219 

-.01637 

.075 

54 

*86399 

-.00603 

• 823 

55 

.90260 

-.00224 

• 7fc7 

56 

•93641 

.00102 

.761 

57 

.96395 

.00198 

.751 

58 

.98400 

.00142 

.757 

59 

.99602 

.00045 

.767 

60 

1,00000 

-•ocooo 

• 773 


AlPHAO ■ 3.B5 DEG® EE S C"0 • -.0909 ETA - 1,131 

VARIABLE GEOMETRY AIRFOIL 664 
DELETED POINTS 52 THROUGH 61 
DELETED POINT 7 

DELETEO POINT * 

DELETED POINTS 2 THROUGH 3 

INSERTED POINT ON UPPTR SURFACE AT X/C ■ 1.1000 

INSERTED POINT ON UPPE® SURFACE AT X/C ■ 1*2000 

INSERTED POINT ON LOWFP SURFACE AT X/C ■ .8400 

INSERTED POINT HN LOWER SURFACE AT X/C • .9600 

INSERTED POINT CN LOWER SURFACE AT X/C • X.ObOO 

INSERTED POINT ON LOWER SURFACE AT X/C • 1.2000 


THE ABOVfc ANGLES OF ATTACK RELATIVE TO THE CHORD LINE 


Y/C ■ -.0350 
Y/C ■ -.0900 
Y/C • -.0265 
Y/C • -.0290 
Y/C ■ -.0490 
Y/C ■ -.0900 
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panel METHOD AIRFOIL 664 CA • 1.58921. 6*76770 AL PH A3 -13.10 DEGREES 

AIRFOIL 664 16.63? THICKNESS 0.00? FLAP 0.00 DEGREES DEFLECTION 

0.00 

N X ■' Y VELOCITY DISTRIBUTIONS FOR THE ABOVE ANGLES OF ATTACK RELATIVE TO THE CHORD LINE 


0 

1.20000 

-.09000 

.019 

1 

1.19373 

-.00620 

.920 

2 

1.17500 

-.07496 

1.021 

3 

1 • 14 ?9 1 

-. 05700 

1,137 

4 

1 . 10 C 30 

-.03500 

1.261 

5 

1.07054 

-.07294 

1 .290 

6 

1.03709 

-.01132 

1.294 

7 

1.00000 

0.00000 

1.205 

9 

.96923 

.00381 

1.286 

9 

.92110 

.02193 

1.294 

10 

.85407 

.03927 

1.303 

11 

.01512 

.04942 

1.326 

12 

.77353 

.06020 

1.359 

13 

.73000 

.07122 

1.401 

14 

*68549 

* 0 P 1 97 

1.455 

15 

• 64 C 43 

.09167 

1.516 

16 

• *9497 

.09937 

1.556 

17 

• 54 069 

. 104 e 2 

1.563 

16 

. 50167 

.10040 

1.568 

19 

.45437 

.11029 

1.572 

20 

,40727 

.11060 

1.577 

21 

.36037 

• 1093 P 

1.562 

22 

.31564 

.10670 

1.589 

23 

.27205 

•10262 

1.596 

24 

♦23051 

.09720 

1.605 

25 

.19145 

.09055 

1.614 

26 

.15521 

.08277 

1.623 

27 

. 12216 

.07401 

1.633 

20 

.09258 

.06441 

1.643 

29 

.06674 

.05416 

1.652 

30 

.04407 

.04348 

1.661 

31 

.02714 

.03261 

1.663 

32 

.01371 

.02183 

1.661 

33 

.00466 

.01155 

1.612 

34 

.00023 

.00229 

1.553 

35 

.00146 

-.00521 

.860 

36 

.00903 

-.01173 

.072 

37 

.02234 

-.01817 

• 263 

30 

• 040^7 

-.02423 

.479 

39 

•06471 

-.02979 

.597 

40 

.09334 

-.03402 

.671 

41 

.12651 

-.03936 

.722 
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AlftFrm 16.631 THICKNESS O.OOt FLAP 0*00 DEGREES Dt FL EC T 1 ON 

0 . 00 

N x t VELOCITY DISTRIBUTIONS FOR THE ABOVE ANGLES OF ATTACK RELATIVE TO THE CHORD LINE 


*2 

.18330 

-.0*3*1 

.759 

A3 

• 20*7* 

- . 0*693 

.787 

A A 

#2*682 

-.0*990 

. 8 08 

A5 

.2955? 

-.05229 

• 62* 

A6 

#3**29 

-.05*06 

.e35 

A7 

.39*5? 

-.05522 

.8** 

A8 

• **556 

-.05572 

.050 

A9 

. *967* 

-.055*8 

.055 

50 

.5*75* 

-.05*33 

*056 

51 

•59719 

-.05219 

.05* 

5? 

.6*51? 

-.0*667 

,0*6 

53 

.69117 

-.0*322 

• 601 

5 A 

.7*1** 

-.036*1 

.7*7 

55 

.79152 

-.030*1 

.700 

56 

.8*000 

-.02650 

.666 

57 

•88180 

-.025*2 

.653 

58 

.92188 

-.02639 

.6*8 

59 

.96000 

-.02900 

.6*7 

60 

1.00*11 

-.03*01 

.6** 

61 

1. 0**21 

-.0*075 

,6*7 

6? 

1 .08000 

-.0*900 

.660 

63 

1.12280 

-.06177 

.698 

6A 

1.15631 

-.07351 

.736 

65 

1.180*7 

-.0825* 

.757 

68 

1.19510 

-.08812 

.77* 

67 

1,20000 

-.09000 

,819 
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BOUNDARY-LAYER DISPLACEMENT ITERATION 

The theoretical results for c x versus from the previous version of 
the computer program (ref. 1) agree remarkably well with experimental measure- 
ments. (For example, see ref. 3.) This good agreement, however, does not hold 
for c x versus a or c m versus a, particularly for aft-loaded airfoils. 
This Is not surprising In that the boundary-! ayer displacement effect was only 
accounted for by reducing the lift-curve slope from Its theoretical value to 
2*. An Improvement could therefore be expected from a more detailed analysis 
of the displacement effect. 

There exists, however, a fundamental flaw In the philosophy of the 
application of displacement Iterations. All displacement effects are of 
second order In boundary-layer theory (ref. 4). Accordingly, It Is Inconsis- 
tent to Include the displacement effect while neglecting other pertinent 
second- order effects which arise from the pressure gradient normal to the 
surface within the boundary layer and other y-component terms In the Navler- 
Stokes equations. This flaw becomes more significant as the boundary-layer 
thickness Increases. 

At the trailing edge, difficult problems arise. The potential-flow 
solution yields steep pressure gradients toward the trailing edqe which 
result In a very high slope for the displacement thickness. This high slope 
can result In a rapid divergence for the displacement Iteration, even for high 
Reynolds numbers. The order (quality) of the tral ling-edge treatment has a 
significant Influence on the results. The wake solution Incorporated In the 
present panel method gives very precise results for the lift coefficients of 
airfoils with blunt trailing edges. It, however, also predicts steep pressure 
gradients toward the trailing edge which. In turn, accelerate the divergence 
of the displacement Iteration. Moreover, this solution also clearly shows 
that the small region which surrounds the trailing edge has a great Influence 
on the solution for the entire airfoil. 

One solution to this divergence problem Is to artificially smooth the 
boundary- layer displacement after each Iteration. But, even If convergence Is 
obtained and, furthermore, even If smoothing were not required for 
convergence, the iteration process would still be questionable due to the 
neglect of the second-order boundary- layer terms previously mentioned. A 
wake solution which minimizes the pressure gradients near the trailing edge 
could Improve the Iteration process but woufd not eliminate the fundamental 
flaw In philosophy. 


10 



The question remains as to what simple procedures can be developed to 
obtain at least a rough estimate of the displacement effect. As previously 
explained, multiple iterations are not logical. Accordingly, in the present 
method, only one Iteration Is performed. The displacement thickness is 
smoothed once and then added to the airfoil contour. The lift and pitching- 
moment coefficients are then computed for the new contour and stored. Later 
the linear portions of the c t - o and % - a curves are adjusted by a 
least-squares fit to these stored values. The separation corrections are then 
applied as discussed in reference 1. Thus, only a few angles of attack require 
this displacement Iteration. The remaining angles of attack are adjusted 
according to the least-squares fit. The displacement effect Is considered to 
be linear In a. A higher-order effect cannot be expected from such a simple 
approach. 


This simple procedure does not require much computing time. The results, 
of course, depend on the smoothing process. In the present version of the pro- 
gram, the curvature of d-j Cx) |l.e.» jls limited. ^|he limit can be 

specified In the input. This limit (SLM) Is preset to i- — jr- < 0.5 (= SLM). 

dx 

The single Iteration Is Initiated by one Input card, which must immedi- 
ately precede an Input card which Initiates a boundary-layer computation 
(l.e., an RE , FLZW, or PLW card). 


DPIT Card 

NUPA, NUPE, NUPI , and NUPU are neglected. 

The F-words specify the angles of attack for which a displacement iter- 
ation Is performed and also the plot mode mbt. The five digits of 
F- are denoted abcde. The variable abc is interpreted as an integer 
n 1 and a displacement Iteration Is Initiated for the nth angle of attack 
on the preceding ALFA card. A displacement Iteration is performed for each 
Reynolds number from the Immediately following RE , FLZW, or PLW card. 

The variable d determines the plot mode. If d > 0, a diagram 
containing the airfoil contour (Including the displacement thickness) and 
the velocity distribution for the angle of attack under consideration Is 
plotted after each displacement Iteration. The plot mode mbt Is set 
equal to d - 1 and Is described under M DIAG Card" In reference 1 (p. 52) 
and reviewed below. 

d = 1 - Axes are drawn, one set of data Is plotted, and the diagram 
Is terminated (l.e., closed to further plotting). 

d * 2 — Axes are drawn, one set of data is plotted, and the diagram 
is open to further plotting. 
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d - 3 — No axes are drawn and one set of data is plotted into the 
existing diagram which is then terminated. 


d * 4 - No axes are drawn and one set of data Is plotted Into the 
existing diagram which remains open to further plotting. 


If d * 2 or 3, the RE , FLZW, or PLW card must specify only one 
Reynolds number. 


Up to five F-words are allowed which means that displacement Iterations 
can be performed for up to five angles of attack. 

d2§i 

If Fc < 0, the limit for — — is set to SLM * -O.bcde. 

5 dx 2 


This new 

limit Is used until It Is reset by another DPIT card with Fc < 0. 
Obviously, only four angles of attack can be specified on DPIT cards 
with F 5 < 0 . 


Examples 

The following card sequences Illustrate some of the OPIT-card options. 


TPftl 0315 1650 400 1750 100 1350 400 3050 4«o 3350 430 3450 470 3650 550 


l*J i 

tpai 
>' > '« 
TRA1 

r* » ~ « 

TRA2 

i» * 

ALT A 


PI AG 


DPIT 

, ~i~ V 

RE 

r» t * 
CDCL 


i i nnii)i)n»iiiM>i<ininnwaaii>»»i)BBx»»ii ■ mm *1 ««« m * m » u u w mm « m m » « om m mv m mn~nlTn w Wim ijj) 

0315 2850 710 3050 1030 HOC 1670 3250 20 3450 70 3*50 90 3850 100 

0315 4050 85 4?50 55 4450 -05 4*50 -195 4350 -2*5 6000 100 

~1 mTm iTimi j ■ II g a 5 S in 5 a i » a aw i h g 5 m ti tf 5 « i « m m w u jj j 5 » 5 »»» i u w ■ g »> j g j n j a » j j n j n ■) 

0315 400 1650 ?00 400 770 600 1 150 200 300 650 300 300 000 000 

iitimiiiMimimwwaiiaini gjfij miniiii»»uiniti<i«iii«MMinynii»m»«Mi)miiminMMiiMi| 

10 000 100 300 300 400 500 600 700 800 900 

i i i mmuimim o»nium)iiMiiatiMiiiMiiiniiiiiiin»iniPMiiiM»iiiminnMii| 

i »i>an mm w » ii S nn aainTn ■■■mbbw h>i«iio<mhimih ww»miwti«owiiwnwMmn»»»| 

110 510 91<I 

i » « an BiimiiiiisiiBai i iiffinyaianiisMiitfflaaiiiaaiitiMtiniiiMisiBniiiiiiiiiwiiHiiMHMs] 

03 200003 6000 

iiidii imia»iiMiimaii»»yiimi«»iii»iiia<n4i)«Minwiiiiiy»MniiMw«MMin«nnB"Biminiii) 
l * n I) n 11 w m II n « i K a a » aim mi)iB>in)>nnima«iin<iinBHMr»M4wBB>MBmiinn)imiiMil 


ENPE 


i i wn iiiaii"iiVuini)inaM»ini«iiaBii»»)iinii4e«M<i4otiuiiiiaiiDMMHBB»iit)'ieaii nnM< „ M ,| 
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After the RE card Is read, displacement iterations for the J^st, fifth, 
and ninth angles of attack (i.e., a * 0°, 4 , and 8 re ]*5 ive . t ° the zero- lift 
line) are performed for both R ■ 2 x 10 4 * 6 and R * 6 x lO®. A diagram is 
plotted for each displacement Iteration. A potential -flow diagram (no dis- 
placement iteration) is also plotted (DIAG card). Thus, one diagram 
10 velocity distributions (DIAG card) (fig. 2) and six diagrams containing one 
velocity distribution each (DFIT card) (fig. 3) are plotted. T J e c, - « 
and c m - a portions of the boundary-layer summary and Its plot (CDCL card) 
are adjusted according to the computed displacement effect. 

It should be noted that each displacement Iteration requires a solution 
from the panel method. Thus, for an airfoil having 61 points, each ^place- 
ment iteration requires approximately 8 seconds CPU time on a Control Data 6600 
computer. 


0315 1650 40 0 1750 100 1 350 400 3050 410 3350 430 8450 470 8650 550 

f, i ~ * t > i i i a n ii ii www i) naan nowa i~» m u »uwni*»w 4i <» o « a ■ « ■ ff - * . * * * n 9 

Tpi^i o?15 2950 710 3050 1030 000 1670 3350 30 3450 70 3650 90 3950 100 

f Ti ■ « i t i » m. 

Tpaj 0315 4050 35 4350 55 4450 -05 4650 -125 4950 -365 6000 100 

TTTTTTJu 

TRA3 0315 400 1*50 200 400 770 600 1 150 200 300 650 300 300 000 000 

rr mimn it h m n » i> m » a n’a a n»ii»»»ii»iiii»m»M<<» *_* ** h miwmim***" ■ jj] 

AL FA 9 300 700 ' 

~ i i > » m ii ii u' w n > ii < a > n B ~HTTrinn>»ii aa»»iiii»«»inM«i ««»n « u m » BeaMmnumnn ■ 

DIAG 1 ___ 

r^TT T > t »n » ayiT'*)imn»vnnn nuvmmu»innmn»mmmtt9mmmmii^mmvmuu»mu mmmumm»mm9mmnnnnn*nn9m^ 

ALFA 10 000 100 200 300 400 500 600 700 300 900 


-T ■ I I » t I »H 

OP IT 440 930 


RE 03 6000 


J7 i i 1 j j j' 1 <n 11 immixoiina n nTh 

CDCL 


4 t 1 1 1 mu unMMiiimni 




e«df 


jT~ nrrrvn umi» ■»■!!■»>■ »i*Tiii» » m» mm ■»*■■«! ■■■»■» ■ ■»»■■■) 

The preceding card sequence plots one diagram (fig. 4) which contains both 
potential -flow and displacement-iteration shapes and velocity distributions 
for a * 3° and 7°. Note that only one Reynolds number Is considered and that 
displacement Iterations are only performed for a = 3° and 7 . The boundary- 
layer sunmary which follows contains the adjustments due to the computed dis- 
placement effect. AC Is the adjusted angle of attack (relative to the zero- 
lift line). 
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SUI*l»«r *I»Fnil 31 5 ANGIE OF ATTACK RELATIVE TO THE 2ERO-LIFT LINE ALPHAO • 2.95 DEGREES 

* INDICATES VFJ.OCITY REDUCTION WITHIN BUBBLE BELOW ,9A 

R ■ 6000000 HU ■ 3 


ALPHA * 

0.00 DEGREES 



1 S TURS S SEP 

CD 

UPPER 

.5749 0,0000 

.0030 

LOWER 

.4107 0.0000 

.0017 

TOTAL 

CL » 0.000 CO • 

.0047 


CM « -.0602 AC » 

.Id 


ALPHA • 

1.00 DEGREES 



1 S TURB S SEP 

CD 

UPPER 

*5676 0.0000 

.0032 

LOWER 

.4066 0.0000 

• 0016 

total 

CL • .110 CO » 

• 0048 


CM • -.0662 AC • 

1.11 

ALPHA • 

2.00 OEGREES 



1 S TUPB S SEP 

CD 

UPPER 

.5980 0.0000 

.0034 

LOWER 

.4069 0.0000 

,0015 

TOTAL 

CL • .220 CD ■ 

.0049 


CM » -.0702 AC • 

2.05 

ALPHA ■ 

3,00 DEGREES 



1 S TUR8 S SEP 

CD 

UPPER 

• 4067 0.0000 

.0036 

LOWER 

.4052 0.0000 

.0015 

TOTAL 

CL • .330 CO * 

• 0050 


CM » -.0722 AC ■ 

2.96 

ALPHA • 

4.00 DEGREES 



1 S TURB S SEP 

CD 

UPPER 

.6147 0.0000 

• 0038 

LOWER 

.4036 0.0000 

• 0014 

TOTAL 

CL • .440 CO • 

.0052 


CM • -.0742 AC - 

3.91 

ALPHA < 

> 5.00 DEGREES 



1 S TUPP S SEP 

CD 

UPPER 

•6569 0.0000 

• 0042 

LOWER 

,4022 0,0000 

• 0013 

TOTAL 

CL • .550 CO • 

• 0056 


CM • -.0762 AC ■ 

4.85 


If 



ALPHA 0 • 2.95 DECKELS 


SUHHAPY AIRFOIL 315 ANGLE OF ATTACK RELATIVE TO THE ZlRO-LIFT LINE 

• I NO 1C AT F S VELOCITY REDUCTION WITHIN BUBBLE 8ELOW .9 A 

R • 6000000 HU • 3 


ALPHA • 

6.00 DEGREES 



1 S TURR S SEP 

CD 

UPPER 

.7«2 3 0*0000 

• 0051 

LOWER 

• 4007 0*0000 

.0013 

TOTAL 

CL • .660 CP « 

.0064 


CH ■ -.0782 AC ■ 

5.78 

ALPHA ■ 

7.00 DEGREES 



i s me s sep 

CD 

UPPER 

.8430 *0006 

.0063 

LOWER 

.3092 0.0000 

• 0012 

TOTAL 

CL • .769 CD • 

.0075 


CH • -.0801 AC • 

6.71 

ALPHA » 

8.00 DEGREES 



1 S TUP 8 S SEP 

CD 

U»PFR 

.9073 .0036 

*0074 

LOWER 

.3970 0.0000 

.0012 

TOTAL 

CL ■ .876 CO ■ 

.0006 


CH • -.0813 AC • 

7.64 

ALPHA • 

9.00 DEGREES 



1 S TURB S SEP 

CD 

UPPER 

.9667 .0071 

.0085 

LOWER 

.3962 0.0000 

.0011 

TOTAL 

Cl • .902 CD • 

.0096 


CH ■ -.0823 AC ■ 

0.38 


xT 



SINGLE ROUGHNESS ELEMENTS 


Recent flight and wind-tunnel experiments indicate that single roughness 
elements such as flap and aileron hinges and poorly faired spoilers signifi- 
cantly degrade the overall performance of an airplane (ref. 5). With the 
previous version of the program (ref. 1 ), only the effect of roughness on 
boundary-layer transition could be considered. Fixed transition points could 
be specified using transition mode 1 or 2 , whereas premature transition due to 
distributed rouqhness or free-stream turbulence could be analyzed using 
transition modes greater than 3. (See "RE Card," ref. 1, p. 56.) 

In the present version of the program, an option has been added which 
allows the analysis of the effect of single roughness elements on a turbulent 
as well as a laminar boundary layer. The method is described in detail in 
reference 5 and reviewed below. 


The increase a «2 of the boundary-layer momentum thickness 62 due t0 a 
single roughness element of height h Is assumed to depend only on the local 

u h h 

roughness Reynolds number R h = — where u h is the x-component of the 


velocity in the turbulent boundary layer at a distance h from the surface. 
For a turbulent boundary layer, the increase of 62 due to the roughness 
element is assumed to be 


A6 


2 


c 


0.15 


u h h 

tnr c 


where c is the airfoil chord and U*, is the free-stream velocity. An expres- 
sion for the velocity u^ is taken from reference 6 and transformed to the 
variables available in tne boundary-layer method. This yields 


^ = \lc~~ 

U > L f 


2.17 In^^-Rt U 6.5 


where U is the local potential -flow velocity, Cf 

/ u « c \ 

friction coefficient, and R| = ~~J is the Reynolds number based on free- 

stream conditions and airfoil chord. In the skin-friction coefficient, t q is 
the shear stress at the wall and p is the air density. 


jis the local skin- 
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If the boundary layer Is laminar at the position of the roughness element, 
transition Is assumed to occur at that position. This Is specified as h * 0 
which acts as a "latest" transition point. Upstream of that position, any 
transition mode except 1 or 2 (fixed transition) is allowed. This approach is 
more logical for many analyses than fixed transition. In front of which no 
other transition criterion Is applied except transition following laminar 
separation. Fixed transition (mode 1 or 2) alone could result in delayed 
transition at some (high) angles of attack - an effect which Is obviously not 
Intended. 

RE Card 

F-words 11-14 contain the data for single roughness elements. These 
words previously only contained the transition points for transition 
modes 1 and 2 (fixed transition). 

If F 14 < 0, F-words 11-14 specify single roughness elements and, 
therefore, transition modes 1 and 2 cannot be used. The five 
digits of Fii - F 14 are denoted abbcc. For F n - F 13 , a Is 
either a blank or 6? For F ]4 , a Is a minus sign (-). The digits 
bb specify the location of tne roughness element x R In percent 
chord. The digits cc which are read as O.cc specify the roughness 
height h in percent chord. Thus, roughness heights can be specified 
over the range 0.0001 <_ h/c <_ 0.0099. F^ and F^ specify roughness 
elements on the upper surface whereas F| 3 and F.i are for the 
lower surface. If x R - 0 Is specified, no roughASss element is 
introduced for that F-word. Thus, 0, 1, or 2 roughness elements can 
be specified on each surface. 

Fn - F ia are read from each RE card which specifies at least one 
Reynolds num&ir. The roughness elements remain in effect until an RE card 
with F 2 t 0 is read. 

Roughness elements can only be analyzed at positions which are actual air- 
foil coordinates. If x R Is specified at an x/c which does not correspond 
to any of the airfoil coordinates, the roughness-element location Is shifted 
to the next airfoil coordinate downstream of x R . If there Is no airfoil 
coordinate close enough to the desired roughness-element location, one can be 
Inserted using a PAN or FXPR card. (See ref. 1.) 

Examp 1 es 

The following RE card specifies two roughness elements on the upper 
surface at x/c * 0.60 and x/c * 0.80, each with a height h/c - 0.0010, and 
one roughness element on the lower surface at x/c * 0.70, with a height 
h/c = 0.0015. 


03 4000 6010 3010 0000-7015 

■mi |, n»»i‘i'»»«»«‘ | ii »«iim"«» | ‘»ii , '»g 
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The following RE card specifies the same roughness elements on the upper 
surface and none on the lower surface. 

RE 03 4000 6010 8010 OOOQ-OUU1 
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Figure 3. - Continued. 








